Appearance of Flat Bands and Edge States in Boron-Carbon-Nitride Nanoribbons 
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Presence of flat bands and edge states at the Fermi level in graphene nanoribbons with zigzag 
edges is one of the most interesting and attracting properties of nanocarbon materials but it is 
believed that they are quite fragile states and disappear when B and N atoms are doped at around 
the edges. In this paper, we theoretically investigate electronic and magnetic properties of boron- 
carbon-nitride (BCN) nanoribbons with zigzag edges where the outermost C atoms on the edges 
are alternately replaced with B and N atoms using the first principles calculations. We show that 
BCN nanoribbons have the flat bands and edge states at the Fermi level in both H2 rich and poor 
environments. The flat bands are similar to those at graphene nanoribbons with zigzag edges, but 
the distributions of charge and spin densities are different between them. A tight binding model and 
the Hubbard model analysis show that the difference in the distribution of charge and spin densities 
is caused by the different site energies of B and N atoms compared with C atoms. 

PACS numbers: 73.20.At, 73.21. Cd, 73.22.-f 



Graphene nanoribbons (GNRs) are nanometer width 
graphene (Gr) sheets [l], and fabricated by cutting 
Gr using electron-beam lithography Q , unzipping carbon 
nanotubes [4|,|5| and synthesized using bottom up process 
0| . The electronic properties of GNRs strongly depend 
on the edge structures. At the zigzag edges, presence 
of the so-called edge states and flat bands at the Fermi 
level were predicted [l|, Q . Recently, the presence of edge 
states in unzipped GNRs were experimentally confirmed 
by Tao et al. [7j. 

In the honeycomb lattice, there are two nonequivalent 
sites in the unit cell, A and B sublattices, and such sub- 
lattice structure plays decisive role for the electronic and 
magnetic properties of the flat bands and edge states. In 
the tight binding model description, the electronic den- 
sity are accumulated on the one sublattice including out- 
ermost C site but identically zero on the other sublattice 
[l|. In the presence of the spin degree of freedom, fer- 
rimagnetic ordering at each site was predicted using the 
Hubbard model calculation and the first-principles cal- 
culations []], H, @|. Therefore, the edge states reflect the 
sublattice structure in the spacial distribution of elec- 
tronic and spin densities. 

On the other hand, boron-carbon-nitride (BCN) 
nanoribbons should show tunable electronic and mag- 
netic properties by B and N doping, because B and N 
atoms act as acceptors and donors in Gr, respectively. 
The electronic and magnetic properties of BCN nanorib- 
bon in which C atoms are partially substituted with BN 
has been investigated by several authors To|- TB| . The 
most of the work handled the hexagonal boron-nitride 
(BN) and Gr hybridized structures and half metallic or 
ferrimagnetic properties along the edges were reported by 
the several authors [Tol-H^. However, they could not ob- 
tain the flat bands and edge states in BCN nanoribbons 



where B and N atoms are doped near the edges because 
such doping disturbs the electronic structures. 

In this letter, we investigate the electronic and mag- 
netic properties of BCN nanoribbons with zigzag edges 
where the outermost C atoms on the edge are alternately 
replaced with B and N atoms as shown in Fig. [T] (a) us- 
ing the first principles calculations. We show that BCN 
nanoribbons in both H2 rich or poor environment have 
the partially flat bands at the Fermi level. We discuss the 
properties of edge states using the tight binding model 
and the Hubbard model calculations. 

We performed the first-principles calculations based 
on the density functional theories within projector aug- 
mented wave method [16l| a nd the local (spin) density 
approximation [L(S)DA] [17j implemented on vasp code 
[18l . LL9j . We considered four different structures as shown 
in Fig. [1] (a) in order to investigate the stability of H ter- 
mination. The cell size in the x-direction was chosen to 
be 2a with a = 2. 446 A which is the lattice constant of 
free standing Gr optimized in LDA. We imposed the vac- 
uum region in the y and z direction about 12 A thick. The 
cutoff energy of plane wave was chosen to be 400 eV and 
the /c-points sampling was chosen to be 12 x 1 x 1 with 
Monkhost-Pack mesh. The geometry was optimized in 
the plane until maximum force becomes less than 10 -3 
eV/A since the out of plane displacements were relaxed 
into the flat structures. 

First, we shall consider the stability of H termination 
in BCN nanoribbons. The formation energy, Ef orm , of 
BCN nanoribbon is defined as 

N c Nu 

E{ OT m = ^BCN ~ 2£?BN 2~^ Gr 2~^ H2 ' ^ 

where £^bcn, ^bn, Eq t and En 2 are the total energies 
of BCN nanoribbon, BN sheet, free standing Gr and 
isolated H2 molecules, respectively. Here, Ac and A^h 
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FIG. 1. (color online) (a) Schematic illustration of the structure of zigzag BCN nanoribbons. In this figure, the red, blue 
and pink circles indicate B, N and H atoms, respectively, (b) The Gibbs free energy measured from that of B,N structure, 
(c) The isosurface plots of spin distribution of BCN nanoribbons. The isovaules are chosen as ±1.0 x 10 _3 /bohr 3 . The black 
circles represent C atoms, (d) The band structures of nonmagnetic BCN nanoribbons within LDA. (e) The band structures of 
magnetized BCN nanoribbons within LSD A. 



are the number of C and H atoms in unit cell, respec- 
tively. Corresponding spin distributions BCN nanorib- 
bons within LSDA are presented in Fig. [T] (c) by isosur- 
face plots with isovalue ±1.0 x 10 _3 /bohr 3 . The same 
plot for GNR with H termination is shown in Fig. [1] (c-v) 
as the reference. The spin distribution at BH,NH, BH,N 
and B,NH structure edges are ferromagnetic while those 
at B,N and GNR edges are ferrimagnetic. Calculated 
inform for several structures with different magnetic con- 
figurations are summarized in Table 1. We found that 
the magnetized BH,NH structure is stable. 

In this comparison, however, effect of environment, i.e., 
temperature and pressure of H2 gas is absent jio| . In or- 
der to clarify the effect, we shall consider the Gibbs free 
energy which is defined as G = ^f orm — A^hMh 2 /2 with 
chemical potential of hydrogen gas fin 2 [20]. In Fig.[T](b), 
calculated Gibbs free energy measured from that of non- 
magnetic B,N structure, AG, is presented as a function 
of /in 2 • We used the stablest state for each structure for 
the calculation of G. We found that magnetized BH,NH 
nanoribbon is stable in H2 rich environment and nonmag- 
netic B,N nanoribbons is stable in H2 poor environment. 

Next, we shall consider the band structures of BCN 



nanoribbons. The band structures of BCN nanoribbons 
within LDA are shown in Fig. Q](d). In Fig. [1] (d-v), the 
band structure of GNRs in the doubled unit cell within 
LDA are shown as the references. We observed the par- 
tially flat bands at the Fermi level around the Y point in 
BCN nanoribbons with BH,NH and B,N structures, while 
the flat bands vanish if B and N atoms are uniformly sub- 
stituted for outermost C atoms [15] . The length of the 
flat bands of BCN nanoribbons in the wavevector space 
is same as that of GNR. For GNR without H termina- 
tion, the dangling bond states of C atoms are partially 
occupied, but fiat bands are open when we do not pas- 
sivate the dangling bonds 21]. On the other hand, for 
BCN nanoribbon with B,N structure, the dangling bond 
states of N is 2.0 eV below, while those of B is 0.7 eV 
above the Fermi level, resulting in the appearance of the 



TABLE I. Calculated results of formation energies at zero 
temperature E{ orra for several magnetic structures. 



£form (eV) 


BH,NH 


BH,N 


B,NH 


B,N 


LDA 


5.56 


9.20 


11.40 


13.29 


LSDA 


5.54 


9.11 


10.79 


13.36 
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fiat bands at the Fermi level. For BCN nanoribbons with 
BH,N and B,NH structures, we could not observe the flat 
bands at the Fermi level but we found doubly degenerate 
bands consisting from 7r-orbital in \k\ < 7r/4a around the 
Fermi level. 

Figure [1] (e) shows the band structures of BCN 
nanoribbons within LSDA. In Fig. [1] (e-v), the band 
structure of GNRs in the doubled unit cell within LSDA 
are also presented as the references. We found that BCN 
nanoribbons are semiconducting except the BH,N struc- 
ture. The band gaps of BCN nanoribbons are 0.155 eV 
for BH,NH structures and 0.189 eV for B,N structures, 
which are smaller than that of GNR 0.201 eV. Although 
the band structures of BCN nanoribbons are similar to 
those of GNRs, the electronic properties and magnetic 
properties of these nanoribbons are quite different as dis- 
cussed below. 

In the remaining part of the letter, we shall analyze 
the flat band and edge states at the Fermi level observed 
in the BH,NH and B,N structures. First, we shall use the 
tight binding model of 7r-electrons and the Hamiltonian 
is given by [22j 

where Ei is an energy of 7r electron at the site z, cj and Ci 
are the creation and annihilation operators of electrons at 
the lattice site z, respectively, stands for summation 
over the adjacent atoms and Uj is the hopping integral 
of 7r electrons from jth atom to ith atom. E{ are Eb, 
Eq and £^n, the site energies at the B, C and N atoms, 
respectively. Following to the Yoshioka et al, we shall 
assume that the hopping integrals are constant regardless 
of the atoms, i.e., tij = t and Eb > 0, Eq — 0, and 
.En = —Eb- [22] 

Calculated band structure of BCN nanoribbon for 
N = 10 and E^/t = 1 is shown in Fig. [2] (a). We repro- 
duced partially flat bands at E = around k = which 
are quite similar to those within the first-principles cal- 
culations. In Figs. 121(b), the charge distributions of BCN 
nanoribbon in the lowest conduction band at k = (i), 
k = 0.27r/2a (ii) and k = 0.47r/2a (hi) are indicated by 
circles of which areas are proportional to the electronic 
charge. At k = 0, most of the electronic charges are ac- 
cumulated at B and N atoms, and remaining electronic 
charges are distributed on sublattices belonging to the 
second outermost C atoms, which are resemble to those 
of the Klein edge [23|. With increase of the wavevec- 
tor k as shown in Figs, (b-ii) and (b-iii), however, the 
electronic charges are distributed over both sublattices. 
Therefore, the edge states at the BCN nanoribbons differ 
from conventional edge states. 

The edge states can be analytically solved as follow: 
We shall consider semi-infinite graphene with zigzag edge 
where outermost C atoms are alternately replaced with 




FIG. 2. (color online) (a) The band structure of BCN 
nanoribbon for N = 10 and E&/t — 1 within the tight bind- 
ing model, (b) The electronic charge distributions in BCN 
nanoribbon for the states indicated in Fig. [2] (a), (c) Schemat- 
ics of semi-infinite Gr sheet with zigzag edge where the out- 
ermost C atoms are alternately substituted with B and N 
atoms. 

B and N as shown in Fig. [2] (c). Let Ai = —77 + 77" 1 and 
A2 = — 1/(77 + 7 ? -1 ) with 77 = e 2%ka . Then, the amplitudes 
of E = state in mth low can be written as ipm, a = 
A^b, = -(£ b /W + Vb, VV 7 = -Ar^B,'and 

ipm,5 = — (-E , b/£)A™ +1, 0b> where is the amplitude at 
B atom. The edge state is defined as an evanescent state. 
Since |Ai| < 1 for \k\ < ir/3a and | A2 1 < 1 in whole 
Brillouin zone, the flat bands exits < k < 7r/3a. As 
shown in Figs. [1] (c-i) , [1] (c-iv), and [2] (a), this result is 
consistent with numerical results. 

Next, we shall discuss the change in the magnetic 
structures in BCN nanoribbons using the Hubbard model 
which is constructed by adding the on-site Coulomb in- 
teraction term to the tight binding model. The Hamilto- 
nian of the system is given by 

U = J2 EiC U Ci ^ ~ 1 Yl c U c 3,v + U J2 ni ^ Ui ^ (3) 

i,cr (i,j),o- i 

where U represents strength of on-site Coulomb inter- 
action. In this study, we shall assume that the on- 
site Coulomb interaction is independent on the chemical 
species of the site i. We adopted the mean-field approx- 
imation of the Hamiltonian Eq. (|3]). 

The band structures for E&/t = 1 are presented in Fig. 
[3] (a) for U/t = 0.5 (i) and for 1.5 (ii). We obtained spin 
split energy bands which agree well those within LSDA 
calculations shown in Figs. [1] (e-i) and (e-vi). In Fig. [3] 
(b) the distributions of spin densities for E&/t = 1 are 
presented by the circles whose areas are proportional to 
the magnitude of spin density for U/t = 0.5 (i) and for 
1.5 (ii). We found that ferrimagnetic order is realized 
for U/t = 0.5, but ferromagnetic order is realized for 
U/t = 1.5. In Fig. [3] (c) Eb and U dependences of spin 
at several sites are presented in color plot. In this fig- 
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FIG. 3. (color online) (a) The band structure of BCN 
nanoribbons with N = 10 using the Hubbard model for 
U/t = 0.5 (i) and for 1.5 (ii). (b) The spin distribution 
for U/t = 0.5 (i) and for 1.5 (ii). In these figures, up and 
down spin density distributions are represented by the yel- 
low (positive) and light blue (negative) circles whose areas 
are proportional to the magnitude of spin density at the site, 
(c) Color plots of Eb and U dependence of s z at several sites 
shown in the bottom of (b-ii). The site indexes are consistent 
with those defined in Fig. [2] (c). In each panel, the solid line 
represents s z — 0. 



ure, s z = is indicated by the white solid line and the 
contours are indicated by the dashed lines. It should be 
noted that spin density at m, a and at m, j3 site equals 
to those at m, 7 at m, 5 site, respectively. From these 
figures, we found that the magnetic order depends on 
U/Eb, i.e., ferromagnetic order is realized for U/Eb <C 1 
and ferrimagnetic order is realized for U/Eb ^> 1. 

Above results within LSDA and the Hubbard model 
suggest that Eb increases and .En decrease when out- 
ermost atoms are terminated by H atoms. To under- 
stand this behavior, we consider the shift of average 
potential at core, AV, by hydrogenation, denned as 
AV = (y(BH,NH) - 4 BH ' NH >) - (V ( B < N ) - 4 B ' N) ), where 
V's are the averaged potentials at the atomic cores and 
Ep's are the Fermi energies. We obtained AV within 
LDA 2.05 eV for B and -1.99 eV for N. As the reference 
we also obtained — 0.09 eV for the outermost C atoms in 
GNR. We found that the core levels shift upward for B 
atoms and downward for N atoms by passivation with H 
atom, resulting in the considerable increase in ^b, since 
the hopping integral t in Gr is 2.7 eV. These large shift of 
the energies 7r-electron in B and N atoms causes change of 



magnetic structure from ferrimagnetic to ferromagnetic 
order in the edge states. 

In summary, we investigated the electronic and mag- 
netic properties of boron-carbon-nitride (BCN) nanorib- 
bons with zigzag edges where the outermost C atoms on 
the edges are alternately replaced with B and N atoms 
using the first principles calculations. We showed that 
BCN nanoribbons in both H2 rich and poor environment 
have the partially fiat bands at the Fermi level around the 
T point. The band structures of BCN nanoribbons are 
quite similar to those in graphene nanoribbons in doubled 
unit cell but the distribution of charge and spin densities 
are different from those in graphene nanoribbons. We 
analyzed the edge states corresponding to the fiat bands 
using the tight binding model and the Hubbard model 
in order to clarify the origin of difference from those in 
graphene nanoribbons. We showed that the origin of such 
different charge and spin density distributions is the dif- 
ferent site energies of B and N atoms compared with C 
atoms. 
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